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ABSTRACT 

We present multi-colour photometry of four field ultracool dwarfs with the triple-beam 
photometer ULTRACAM. Data were obtained simultaneously in the Sloan-</ band and 
a specially designed narrow-band Nal filter. The previously reported 1.8-hr period of 
Kelu-1 is here recovered in the g'-band, but the lack of any significant variability in 
the Nal light of this object precludes any conclusion as to the cause of the variability. 
2MASS 2057-0252 and DENIS 1441-0945 show no convincing evidence for variability. 
2MASS 1300+1912, on the other hand, shows good evidence for gradual trends in both 
bands at the 5% level. These trends are anti-correlated at a high level of significance, 
a result which is incompatible with models of starspot-induced variability. It would 
seem likely that dust cloud "weather" is responsible for the short-term variability in 
this object. 
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1 INTRODUCTION 

The photometric variability of brown dwarfs and very low- 
mass stars (collectively known as ultracool dwarfs or UCDs) 
has received substantia l attention over the last five years. 
The recent review by Bai ler- Jonesl J2004) suggests that 
around 40% of UCDs which are surveyed show variability. 
The variability appears to be transient in nature - for ex- 
ample, obser v ations of t he L3 dwarf 2MASS 1146+2230 by 
iGelino et all J2002D and IClarke et all J2002D failed to fin d 
the variability reported by TB^tile^^ones fc Mundtl i200ll) . 
Also, e xtensive observations of the M9.5 dwarf BRI 0021- 
0214 bv lMartm et alJ (1200 lT) found both periodic and tran- 
sient /-band variability that did not correspond to the ex- 
pected rotation period. It is therefore possible that aZZUCDs 
show variability at some level. 

The variability falls into two broad categories: periodic 
modulations with periods up to a few days and magni- 
tudes of up to 100 mmag, and non-periodic variations with 

1 Based on observations made at the European Southern Obser- 
vatory, Paranal, Chile (ESQ program 075.C-0505) 



timescales on the order of hours to days and amplitudes of 
~ 10 to ~ 100 mmag in the /-band, where most of the ob- 
servations to date have been undertaken. The existence of 
non-periodic variability is interesting, as the monitoring sur- 
veys undertaken would have been sensitive to the rotation 
per iods of UCDs. The existence of non-periodic variability 
led iBailer- Jones fc Mundtl i200ll) to propose that the vari- 
ability is due to surface features which evolve on timescales 
shorter than the rotational period. 

There are two prime candidates for the surface features 
responsible for variability. Cool, magnetic star-spots are an 
attractive explanation for the very young UCDs observed 
in clusters. Star-spots are known to cause variabilit y in the 
more massive T Tauri stars (see lHerbst et al1ll994l for ex- 
ample), and furthermore, young brown dwarfs show simi- 
lar X-ray properties to T- Tauri stars, implying that a simi- 
lar magnetic activity mechanism is at work in young UCDs 
and T Tauri stars (e.g. lOzawa et al]l2005t lBouvll2004l) . The 
older, field, UCDs may not show magnetic spots, as the 
cooler photospheres will be neutral and re sult in a weak 
coup l ing between the gas a nd magnetic field f Mohantv et alJ 
l20oilFleming et al-lEoOol) . For the field UCDs, dust clouds 
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Figure 1. Calculations of the change in the spectrum of a UCD 
with T e ff = 1900K, due to the appearance of a cloud or spot 
which covers 10% of the stellar surfa ce. The c alculations were 
made using the synthetic spectra of Alla rd et all bOOll) , The 
four lines represent a the change in the spectrum of the UCD, 
caused by the formation of: a clear hole in a dusty atmosphere 
(dashed line), a dusty cloud on a clear atmosphere (solid line), 
a 200K cooler spot on a dusty atmosphere (intermittent dashed 
line) and a 200K cooler spot on a clear atmosphere (dotted line). 
Also shown are the bandpasses of the Sloan-g' filter, and the spe- 
cially designed Nal filter (A cen = 5900A, <5A = 220A) used for the 
observations presented in this paper. 



are a more plausible candidate: dust should form under 
the cool and dense conditions found in their photospheres, 
and dusty atmospheric models are more successful in repro- 
ducing the observed colors of UCDs tha n dust-free models 
^Leggett et al.lll998l : iMartfn et all 1200(1 ). The atmospheres 
of UCDs are also likely to be highl y dynamic, as a re sult 
of turbulence excited by convection IjAllard et al]l200ll) . It 
is quite likely that rapidly evolving dust cloud "weather" is 
responsible for the non-periodic variability seen in UCDs. 

The different spectral signatures of variability caused 
by spots and dust clouds allows one, in principle, to dis- 
tinguish between these two candidates based upon time- 
resolved spectrophotometry. Unfortunately, the results so 
far have been inconclusive, as no significant variabili ty was 
detected in dust-sensiti ve wavelength regions dClarke et alJ 
2003 : iBailer- Jonesir2 002 ) . The reasons for this lack of success 
are not clear. It may be simply that variability in UCDs is 
transient, and the limited number of observations to date 
have merely been unlucky. Alternatively, it may be that 
the systematic errors in spectrophotometry, arising from sky 
subtraction, telluric correction and slit-loss correction are 
too large to attain the required accuracies of a percent or 
better implied by the magnitude of variability in the /-band. 

A better approach may be to use time-resolved pho- 
tometry, taken simultaneously in a narrow band centred on 
dust-centred features and a continuum filter. Photometry 
does not suffer from slit-loss errors, as other stars in the 
field can ve used to correct for the light that falls outside 
the aperture. Also, the photons are collected in fewer pixels, 
allowing greater signal to noise ra tios to be achieved. Nar - 
row band photometry was used bv lTinnev fc Tollevl dl999l) . 
who obtained quasi-simultaneous narrow band photometry 
of two UCDs, DENIS-P J1228-1547 and LP 944-20. Vari- 
ability was found in both bands for LP 944-20, showing 



the promise of photometry for future studies. Unfortunately, 
the complex na ture of the spectral beha viour in the narrow 
bands chosen bv lTinnev fc Tollevl ill 9991) . precluded any firm 
conclusion as to the origin of the variabilit y. We outline here 
a mod ification of the method proposed bv lTinnev fc Tollevl 
(1999). Multi-band photometry has also been used succes- 
fully to show that, in some cases at le ast, the varia bility 
of M - dwarfs is due to magnetic spots jRockenfeller et all 
2006). Figure Q shows the d etails of the method (which i s 
described in more detail by Bailer- Jones fc MunddfeoOll) . 
The opacity within a spot feature shows a broadly similar 
colour dependence to the opacity of the immaculate pho- 
tosphere, but the cool spot lowers the average temperature 
of the visible portion of the star. Therefore, the evolution 
or appearance of a spot feature produces a dimming of the 
spectrum at all wavelengths. If we choose to observe simul- 
taneously in two filters, one centred on the continuum, the 
other centred on the neutral alkali lines of Nal, a cool spot 
will cause magnitude changes which are correlated between 
the filters. In the case of the formation of a clearing within a 
dusty cloud deck, the clearing removes opacity in the contin- 
uum, but significantly increases opacity in spectral regions 
dominated by molecular lines (the strong resonance lines of 
Nal seen at 0.59 /im are good examples of such a region). 
A dust cloud forming on a clear photosphere has the oppo- 
site effect, increasing the continuum opacity, but decreasing 
the opacity at Nal wavelengths. Thus, in contrast to the 
case of a starspot, a change in the cloud coverage of a UCD 
causes anti- correlated magnitude changes between the two 
filters. By centering our narrow band filter on the neutral 
alkali lines, we can thus distinguish between spots and dust 
clouds. 

In this paper we present multi-colour photometry of 
four field ul tracool dwarfs with the triple-beam photometer 
ULTRACAM llDhillon fc Marshll200ll) . Data were taken simul- 
taneously in the Sloan-p' band and the specially designed 
narrow-band Nal filter shown in figure The observations 
are described in section |5J the results presented in section [3] 
and discussed in sectional whilst in section |5] we draw our 
conclusions. 



2 OBSERVATIONS 

The UCDs 2MASS 1300+1912, Kelu-1, DENIS 1441-0945 
and 2MASS 2057-0252 were observed simultaneously in the 
Sloan-p' and Nal bands using ULTRACAM on the 8.2-m meli- 
pal unit of the Very Large Telescope (VLT) at Paranal, 
Chile. Thanks to the frame-transfer CCDs employed in UL- 
TRACAM we were able to obtain 10 second exposures with no 
dead-time between frames. Sloan-u' images were obtained 
simultaneously with the other colours, but the UCDs are 
too faint to produce useful lightcurves in this band. The ob- 
servations are summarised in Table Q Data reduction was 
carried out using the ULTRACAM pipeline date reduction soft- 
ware. Extraction of target and comparison star ligh tcurves 
was p erformed using an optimal extraction method iNavlorl 
1998). Due to the 2. 6' field of view of ULTRACAM on the VLT 
and the sparse nature of the target fields, a limited number 
of comparison stars are available for each object. 

On inspection of the images, it is clear that the im- 
ages suffer from structured vignetting and scattered light 
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Table 1. Journal of Observations. 



UT Date Target UT UT Seeing Airmass 

Start Start End arcsecs 



2005-05-08 2MASS 1300+1912 23:22 04:13 0.6-1.0 1.4-2.3 

2005-05-13 DENIS 1441-0945 03:11 07:02 0.5-0.8 1.0-1.4 

2005-05-13 2MASS 2057-0252 07:05 10:14 0.5-0.8 1.1-1.5 

2005-05-14 Kelu-1 02:37 06:49 0.5-1.0 1.0-2.0 

2005-05-17 2MASS 2057-0252 07:52 09:50 0.8-1.1 1.0-1.1 



components at the 10% level. Further investigation of these 
effects suggest that they arise in the collimator optics. To 
prevent the scattered light in the flat-field frames from in- 
troducing an error in the flat-field process, large-scale trends 
were removed from the fiat-field before use, by dividing the 
flat-field by a median-filtered version of itself. This means 
that vignetting is not corrected for by our flat field. For the 
majority of our observations this is not a major problem, 
as the position of the stars on the chip remained stable to 
within 1 pixel throughout observation. Whilst the absolute 
flux level for these stars will be in error because of the lack of 
vignetting correction, the shape of the light-curves will not 
be affected. The one exception to this is the observation of 
2MASS 1300+1912. At various times in the observation of 
this object, the target was moved around the chip. In prin- 
ciple this might introduce spurious features into the light 
curve, as the relative vignetting between target and compar- 
ison star changes. In practice, the distance between compar- 
ison star and target was small (~140 pixels), and we believe 
this effect is negligible, a conclusion which is reinforced by 
the relative flatness of the comparison star lightcurve in fig- 
ure |21 and the lack of sudden jumps in the lightcurves at 
times when the star was moved. 



3 RESULTS 

3.1 Rapid Variability 

We performed a search for rapid variability in all our objects, 
by calculating the reduced x _st l uare d with respect to the 
polynomial fits shown in figure |5] In each case, the reduced 
X-squared with rspect to the fits was less than or equal to 
a value of 1.5. We conclude that no objects show evidence 
for "rapid" variability, i.e variability on timescales of a few 
minutes. In the rest of the paper, the term variability refers 
to medium term trends within the data, i.e variability on 
the timescale of hours. 



3.2 2MASS 2057-0252 

2MASS 2057-02 52 is an LI . 5 dwarf with an /-band mag- 
nitude of 16. 5 llKoenll2003ft . The object showed tentative 
evidence for I-band variability on one nig ht, which w as ab- 
sent on the second night it was observed (lKoenll2003l) . Our 
observations of 2MASS 2057-0252 are presented in figure |5] 
No evidence of variability is present in the Nal lightcurve, 
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Figure 3. The (/-band lightcurve of Kelu-1, phased on the known 
1.8-hr period, and binned into nineteen phase bins. Data taken 
at high airmass and with a poor signal-to-noise ratio (after MJD 
53504.22) was masked out prior to phasing. 



as determined by a one-sided F-test. Whilst the (/-band 
lightcurve shows a steady decline in flux with time, the flux 
is strongly correlated with airmass, which strongly suggests 
that this is due to a colour effect, arising because the target 
is redder than the comparison star. This also naturally ex- 
plains the lack of variability in the Nal filter, as this filter 
is much narrower, and redder, than the Sloan-g' filter. We 
therefore conclude that this object shows no strong evidence 
for variability on a 3 hour timescale. 



3.3 Kelu-1 



Kelu-1 is an L2 dwarf llRuiz et alJll997ft, and a known bi- 
nary star (ILiu fc Leggetj 120051 : iGelino et alJl2005i) . Obser- 
vations centred on the complex of molecular bandheads at 
■-^8600 A show strong per iodic variability, with a period of 
1.8-hr JClarke et alJ l2002D . Time-resolved spectroscopy of 
the object showed no evidence for variability in the dust- 
sensitive TiO bandheads, however JClarke et alj E003'l. Our 
observations are presented in figure |2] The Nal lightcurve 
shows no evidence of variability: we used a one-sided F-test 
to show that a l s *-order polynomial was a better match to 
the data than a constant fit at a significance of just 50%. 

The (/-band flux, however, shows a rising trend with 
time. This rising trend is strongly correlated with airmass 
and is most likely a colour effect. Additional variability is 
visible, however, at the start of the observations, where the 
airmass was low, and signal-to-noise is better. A Lomb- 
Scargle periodogram of the (/-band data taken before MJD 
53504.22 (airmass 1.3) has a peak at the known 1.8-hr period 
of Kel u-1. Following the analysis of ISchwarzenberg-Czernvl 
( 1998), the significance of this peak is better than 99.9 per- 
cent. The formal significance should be treated with caution; 
non- white noise, deviations of the lightcurve from a pure si- 
nusoid and irregular sampling can all affect the periodogram, 
and the length of our run is only just longer than the ex- 
pected period. However, the significance is high enough for 
us to conclu de that the 1. 8 -hr p eriodic variability at 8600A 
reported by IClarke et alJ 12002^1 is also present in the g'- 
band light of Kelu-1. The phased lightcurve is plotted in 
figure [3] 
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Figure 2. Clockwise from top left: Lightcurves of 2MASS 2057-0252, Kelu-1, 2MASS 1300+1912 and DENIS 1441-0945. The flux of 
the target star, divided by a nearby comparison, is plotted against Modified Julian Date. The lightcurves have been binned to a time 
resolution of approximately one minute, normalised by division of the mean value, and a vertical offset of 0.2 has been applied to the 
Sloan-</ lightcurve. The Nal lightcurvc is plotted with red circles, the Sloan-g' lightcurve is plotted with green triangles. The lightcurves 
of the comparison star, divided by a third, nearby comparison star are plotted beneath the target star's lightcurves, with the same 
symbols. These lightcurves have also been normalised and an offset added for clarity. Also shown are I s * order polynomial fits to the 
UCD data (except in the case of DENIS 1441-0945, where a 2 nd order fit is shown). Dashed vertical lines in the 2MASS 1300+1912 
panel show the times at which the star's position on the chip was changed. For each object additional panels show the seeing in arcsecs 
(middle panel) as measured from the FWHM of stellar images, and the airmass of the UCD (top panel). 



3.4 2MASS 1300+1912 

2MASS 1300 +1912 is an LI dwa rf with an /-band magni- 
tude of 15.9 iGelino et al . 2002). Convincing evidence for 
non-periodic J-band variability with an amplitude of ~0.2 



mags was reported bv lGelino et alJ i2002l) . Our observations 
are presented in figure |2] The observations taken above air- 
mass of 1.7 are not shown as the signal/noise ratio is too 
low to be useful. It is clear from figure [5] that the object 
shows variability at a level of about 5 percent in both bands. 
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This variability is significant, and clearly anti-correlated 
(the correlation coefficient between the two lightcurves is 
—0.09 ± 0.01). To determine the significance of the anti- 
correlation we computed 2 nd order polynomial fits to the 
lightcurves of the target and comparison star. The differ- 
ence between the (/-band and Nal slopes is —0.6 ± 0.1 and 
0.00 ± 0.02 for the target star and comparison star respec- 
tively. Futhermore, the individual slopes of the comparison 
star's lightcurves are both consistent with zero. We con- 
clude that the anti-correlation is significant, and not intro- 
duced by variability in the comparison star. Although the g- 
band lightcurves of the other objects are affected by colour- 
term effects, this is not the case here. Whilst the airmass 
first falls and then rises, the gr'-band lightcurve of 2MASS 
1200+1912 shows a steady rising trend. We conclude that 
colour-term effects do not affect the gr'-band lightcurve of 
2MASS 1300+1912. Furthermore, the variability is not cor- 
related with the seeing, sky brightness or transparency, and 
we conclude that systematic errors resulting from the re- 
moval and/or correction of these factors do not cause the 
observed variability. 

Such variability, if real, would be strong evidence that 
the variability in this system is caused by dust clouds. Un- 
fortunately, the position of the target on the CCD was not 
constant throughout the run. As discussed in section |5] this 
could lead to spurious variability in the lightcurve. A rea- 
sonable check of the presence of such effects would be the 
structure in the differential lightcurve of our chosen com- 
parison and another star on the CCD. Only one such star 
is present, approximately 530 pixels away from our compar- 
ison star. The differential lightcurve of our comparison star 
divided by this second star is plotted in figure [5] Variability 
is indeed present, but only at a level of approximately 1%. 
Futhermore, there are no sudden jumps in the lightcurves 
at the times when the star's position was changed. We con- 
clude that vignetting has not introduced features into the 
lightcurve, although the variability in the lightcurve at the 
1% level suggests that one of the comparison stars is itself 
variable, albeit at a low level. 

We therefore conclude that the anti-correlated variabil- 
ity seen in 2MASS 1300+1912 is real. This conclusion is 
strengthened by the fact that most systematic errors would 
produce spurious variability that is correlated between the 
two wavebands. In particular, this applies to both colour 
effects, and any errors introduced by the uncorrected vi- 
gnetting (the vignetting patterns are very similar in both 
CCDs). 



3.5 DENIS 1441-0945 

DENIS 1441-0945 is an LI dwarf with an I- band magnitud e 
of 16.9 (Koen 2003), and a know n binary l|Bouv et all 2003'). 
It was observed in the /-band bv lKoenl (120031) - no convinc- 
ing evidence for variability was reported. Although a rising 
trend was present in the lightcurve, this may have been a 
colour-term effect. Our observations of DENIS1441-0945 are 
presented in figure [5] No evidence of variability is present 
in the Nal lightcurve, as determined by a one-sided F-test. 
Whilst the (/'-band lightcurve shows a gradual change in 
flux with time, the flux is strongly correlated with airmass, 
which suggests that this is due to a colour-term effect. We 



therefore conclude that this object shows no evidence for 
variability. 



4 DISCUSSION 

The detection of anti-correlated variability between the g'- 
band and the Nal band in 2MASS 1300+1912 is incompat- 
ible with models of starspot-induced variability. Given that 
anti-correlated variability is expected from models of dust- 
cloud variability, and that rapid evolution of dust clouds is 
to be expected in UCDs, it is likely that dust clouds are re- 
sponsible for the variability in this system. The relative am- 
plitude of variability between the two bands can constrain 
the nature of the cloud evolution responsible for variability. 
The models used in figure show that a dust cloud forming 
on a predominantly clear amosphere would produce a much 
stronger signal in the Nal band than the Sloan-gr' band. In 
contrast, a small clearing in a predominantly dusty atmo- 
sphere produces variability of about the same magnitude in 
each band. The amplitude of variability can also constrain 
the size of the cloud or clearing; our observations are consis- 
tent with the development of a clear patch covering approx- 
imately 4 percent of a predominantly dusty atmosphere. 

The result presented here demonstrates that simulta- 
neous, multi-colour photometry is effective in distinguishing 
between different variability scenarios for UCDs, and opens 
up the possibility of understanding the variability of UCDs 
in more detail. Clearly, the method should be applied to a 
large sample of UCDs to determine if the variability in all 
UCDs is due to dust clouds, and at what effective tempera- 
tures dust-cloud variability becomes important. 



5 CONCLUSIONS 

We present simultaneous, multi-colour photometry of four 
ultracool dwarfs with the triple-beam photometer ultra- 
CAM. Data were obtained simultaneously in the Sloan-p' 
band and a specially designed narrow-band Nal filter. Of 
the four objects, only 2MASS 1300+1912 shows good ev- 
idence for variability in both bands. For this object, the 
variability is anti-correlated at a high level of significance, 
providing the first direct evidence that dust-cloud weather 
is responsible for the variability in ultracool dwarfs. 
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